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Abstract 

We develop a multi-curve term structure setup in which the modelling ingredients 
are expressed by rational functionals of Markov processes. We calibrate to LIBOR 
swaptions data and show that a rational two-factor lognormal multi-curve model is suf¬ 
ficient to match market data with accuracy. We elucidate the relationship between the 
models developed and calibrated under a risk-neutral measure Q and their consistent 
equivalence class under the real-world probability measure P. The consistent P-pricing 
models are applied to compute the risk exposures which may be required to comply with 
regulatory obligations. In order to compute counterparty-risk valuation adjustments, 
such as CVA, we show how positive default intensity processes with rational form can 
be derived. We flesh out our study by applying the results to a basis swap contract. 

Keywords: Multi-curve interest rate term structure, forward LIBOR process, rational 
asset pricing models, calibration, counterparty-risk, risk management, Markov func¬ 
tionals, basis swap. 


1 Introduction 

In this work we endeavour to develop multi-curve interest rate models which extend to 
counterparty risk models in a consistent fashion. The aim is the pricing and risk manage¬ 
ment of financial instruments with price models capable of discounting at multiple rates 
(e.g. OIS and LIBOR) and which allow for corrections in the asset’s valuation scheme so 
to adjust for counterparty-risk inclusive of credit, debt, and liquidity risk. We thus propose 
factor-models for (i) the Overnight Index Swap (OIS) rate, (ii) the London Interbank Offer 
Rate (LIBOR), and (iii) the default intensities of two counterparties involved in bilateral 
OTC derivative transactions. The three ingredients are characterised by a feature they 
share in common: the rate and intensity models are all rational functions of the underlying 
factor processes. In choosing this class of models, we look at a number of properties we 
would like the models to exhibit. They should be flexible enough to allow for the pricing of 
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a range of financial assets given that all need discounting and no security is insulated from 
counterparty-risk. Since we have in mind the pricing of assets as well as the management 
of risk exposures, we also need to work within a setup that maintains price consistency 
under various probability measures. We will for instance want to price derivatives by mak¬ 
ing use of a risk-neutral measure Q while analysing the statistics of risk exposures under 
the real-world measure P. This point is particularly important when we calibrate the in¬ 
terest rate models to derivatives data, such as (implied) volatilities, and then apply the 
calibrated models to compute counterparty-risk valuation adjustments to comply with reg¬ 
ulatory requirements. The presented rational models allow us to develop a comprehensive 
framework that begins with an OIS model, evolves to an approach for constructing the 
LIBOR process, includes the pricing of fixed-income assets and model calibration, anal¬ 
yses risk exposures, and concludes with a credit risk model that leads to the analysis of 
counterparty-risk valuation adjustments (XVA). 

The issue of how to model multi-curve interest rates and incorporate counterparty-risk 
valuation adjustments in a pricing framework has motivated much research. For instance, 


research on multi-curve interest rate modelling is presented in Ivijima, Tanaka, and Wong 


(2009j), Kenyon (2010), Henrard 

(2007, 

2010, 

2014 

), Bianchetti (2010), Mercurio (2010b] 

2010a, 2010c 

), Fujii, Shimada, and Takahashi 

(2011, 2010), Moreni and Pallavicini (2014), 


Bianchetti and Morini (2013), [Filipovic and Trolle (2013) and Crepey, Grbac, Ngor and 


Skovmand (2014). On counterparty-risk valuation adjustment, we mention two recent books 


by Brigo, Morini, and Pallavicini (2013) and Crepey, Bielecki and Brigo (2014); more refer¬ 
ences are given as we go along. Pricing models with rational form have also appeared before. 


Flesaker and Hughston (1996) pioneered such pricing models and in particular introduced 
the so-called rational log-normal model for discount bond prices. For further contribu¬ 


tions and studies in this context we refer to Rutkowski (1997), Doberlein and Schweizer 


(2001) and Hunt and Kennedy (2004). More recent work on rational pricing models include 


|Brody and Hughston (2004), Hughston and Rafailidis (2005), Brody, Hughston and Mackie 


|(2012|), |Akahori, Hishida, Teichmann and Tsuchiya (2014|), Filipovic, Larsson and Trolle 


|(2014| ), |Macrina and Parbhoo (2014 ), and Nguyen and Seifried (2014). However, as far as 
we know, the present paper is the first to apply such models in a multi-curve setup, along 
with Nguyen and Seifried (2014), who develop a rational multi-curve model based on a 
multiplicative spread. It is the only one to deal with XVA computations. We shall see that, 
despite the simplicity of these models, their performances in these regards are comparable 
to those by Crepey, Grbac, Ngor and Skovmand (2014]) or Moreni and Pallavicini (2013 


2014). Other recent related research includes Filipovic, Larsson and Trolle (2014), for 


the study of unspanned volatility and its regulatory implications, Cuchiero, Keller-Ressel 


and Teichmann (2012), for moment computations in financial applications, and Cheng and 


Tehranchi (2014), motivated by stochastic volatility modelling. 


We give a brief overview of this paper. In Section [2] we introduce the rational models 
for multi-curve term structures whereby we derive the forward LIBOR process by pricing a 
forward rate agreement under the real-world probability measure. In doing so we apply a 
pricing kernel model. The short rate model arising from the pricing kernel process is then 
assumed to be a proxy model for the OIS rate. In view of derivative pricing in subsequent 
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sections, we also derive the multi-curve interest rate models by starting with the risk-neutral 
measure. We call this method “bottom-up risk-neutral approach”. In Section[3j we perform 
the so-called “clean valuation” of swaptions written on LIBOR, and analyse three different 
specifications for the OIS-LIBOR dynamics. We explain the advantages one gains from the 
chosen “codebook” for the LIBOR process, which we model as a rational function where the 
denominator is in fact the stochastic discount factor associated with the utilised probability 
measure. In Section [4j we calibrate the three specified multi-curve models and assess them 
for the quality of fit and on positivity of rates and spread. We conclude by singling out 
a two-factor lognormal OIS-LIBOR model for its satisfactory calibration properties and 
acceptable level of tractability. In Section [5] we price a basis swap in closed form without 
taking into account counterparty-risk, that is we again perform a “clean valuation”. In 
this section we take the opportunity to show the explicit relationship in our setup between 
pricing under an equivalent measure and the real-world measure. We compute the risk 
exposure associated with holding a basis swap and plot the quantiles under both probability 
measures for comparison. As an example, we apply Levy random bridges to describe the 
dynamics of the factor processes under P. This enables us to interpret the re-weighting of 
the risk exposure under P as an effect that could be related to, e.g., “forward guidance” 
provided by a central bank. In the last section, we present default intensity processes with 
rational form and compute XVA, that is, the valuation adjustments due to credit, debt, 
and liquidity risk. 


2 Rational multi-curve term structures 

We model a financial market by a filtered probability space (fi, J 7 , P, {PijtKt), where 
P denotes the real probability measure and {J~t\()<t is the market filtration. The no¬ 
arbitrage pricing formula for a generic (non-dividend-paying) financial asset with price 
process { Str}o<t<T , which is characterised by a cash flow Stt at the fixed date T, is given 
by 

S tT = (2.1) 

7 Tt 

where {TTt}o<t<u is the pricing kernel embodying the inter-temporal discounting and risk- 
adjustments, see e.g. Hunt and Kennedy (2004). Once the model for the pricing kernel 


is specified, the OIS discount bond price process {PtT}o<t<T<u is determined as a special 


case of formula (2.1) by 


P tT = | P t \. 

7 Tt 

The associated OIS short rate of interest is obtained by 


( 2 . 2 ) 


n = — {dr lnPtr) \r=t, 


(2.3) 


where it is assumed that the discount bond system is differentiable in its maturity parameter 
T. The rate {ry} is non-negative if the pricing kernel { 7 r*} is a supermartingale and vice 
versa. We next go on to infer a pricing formula for financial derivatives written on LIBOR. 


In doing so, we also derive a price process (2.6) that we identify as determining the dynamics 
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of the forward LIBOR or, as we shall call it, the LIBOR process. It is this formula for the 
LIBOR process that reveals the nature of the so-called multi-curve term structure whereby 
the OIS rate and the LIBOR rates of different tenors are treated as distinct discount rates. 


2.1 Generic multi-curve interest rate models 


We derive multi-curve pricing models for securities written on the LIBOR by starting with 
the valuation of a forward rate agreement (FRA). We consider 0<f<To<T2<---< 
Ti < • • • < T n , where To, T), ..., T n are fixed dates, and let IV be a notional, K a strike rate 
and 5i = Ti — Tj_The fixed leg of the FRA contract is given by NK5i and the floating 
leg payable in arrear at time T) is modelled by -/V<5jT(T); T)_i, Ti), where the random rate 
L(Ti] Tj_i, Tj) is -measurable. Then we define the net cash flow at the maturity date 
Tj of the FRA contract to be 


by 


H T .=N5 i [K-L(T i -,T i -. 1 ,T i 


H t Ti = — E P [irTiHTi | T t ] 

^Ft 


:)]• 

(2.4) 

(|2.1[), that is, 

for 0 < t < T- 1 , 


(2.5) 

, T) | T t ] • 

(2.6) 


where we define the (forward) LIBOR process by 

L(f;T_i,T) := -E p [n Ti L{Ti]l 

The fair spread of the FRA at time t (the value K at time t such that HtTi = 0) is then 
expressed in terms of L(t: T;-i, T)) by 

T(f;T_i,T) 


K, = 


p m 


(2.7) 


For times up to and including T)_i, our LIBOR process can be written in terms of a 
conditional expectation of an TV;., -measurable random variable. In fact, for t < T)_i, 

E p [vrjvL(Tj; T^Ti) \ T t ] = E p [e p [n Ti L{T f , T^, TJ \ Tt^) \ Ft] (2.8) 


= E 


E p \7r Ti | T(Tj;Tj-i,Ti) | F t , (2-9) 


and thus 


1 


L(t, Tj_!, Tj) = —E p E p [t t Tz | T Ti _ ± ] L(T f} T^T,) \ T t 

TTt L 


( 2 - 10 ) 


The (pre-crisis) classical approach to LIBOR modelling defines the price process {H t T, } 
of a FRA by 

H m = N [(1 + SiK)P tTi - P mi _ 1 ] , (2.11) 


see, e.g., Hunt and Kennedy (2004). By equating with (2.5), we see that the classical 
single-curve LIBOR model is obtained in the special case where 


1 


= — (PtTi _! - Pm ) 


( 2 - 12 ) 
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Remark 2.1. In normal market conditions, one expects the positive-spread relation 
L(t;T,T + 6i) < L(t;T,T + Sj), for tenors Sj > Si, to hold. We will return to this re¬ 
lationship in Section [4] where various model specifications are calibrated and the positivity 
of the spread is checked. LIBOR tenor spreads play a role in the pricing of basis swaps, 
which are contracts that exchange LIBOR with one tenor for LIBOR with another, differ¬ 
ent tenor (see Section [5]). For recent work on multi-curve modelling with focus on spread 
modelling, we refer to |Cu chiero, Fontana and Gnoatto (2014]). 


2.2 Multi-curve models with rational form 


In order to construct explicit LIBOR processes, the pricing kernel and the random 


variable L(R,;: , T t ) need to be specified in the definition (2.6). For reasons that will 


become apparent as we move forward in this paper, we opt to apply the rational pricing 


models proposed in Macrina (2014). These models bestow a rational form on the price 


processes, here intended as a “quotient of summands” (slightly abusing the terminology 
that usually refers to a “quotient of polynomials”). This explains the terminology in this 
paper when referring to the class of multi-curve term structures, or to generic asset price 
models, and later also to the models for counterparty risk valuation adjustments. 

The basic pricing model with rational form for a generic financial asset (for short “ra¬ 
tional pricing model”) that we consider is given by 


c __S 0T + b 2 (T)A ( t 2) + b 3 (T)A ( t 3) 

^tT /-|\ 5 

fbt + &i(*)4 } 


(2.13) 


where Sqt is the value of the asset at t = 0. There may be more M-terms in the numerator, 
but two (at most) will be enough for all our purposes in this work. For 0 < t < T and 
i = 1,2,3, b{(t) are deterministic functions and A^ = Ai{t,x\ 1 ^) are martingale processes, 
not necessarily under P but under an equivalent martingale measure M, which are driven 
Marko v pr ocesses {xj 1 ' 1 }. The details of how the expression (1 


by 


the formula (2.1), and in particular how explicit examples for {A t 


(0 


2.13) is derived from 


can be constructed, 
Macrina (2014). Here we only give the pricing kernel model associated with 


are shown in 


the price process (2.13), that is 


TTi = 


M 0 L 


0t 


+ bi{t)A^ 


M t 


(2.14) 


where {AR} is the P-martingale that induces the change of measure from P to an auxiliary 
measure M under which the {A y } are martingales. The deterministic functions Rot and 
b\{t) are defined such that Rot + b\(t)Af 1 ' > is a non-negative M-supermartingale (see e.g. 


Example 2.1), and thus in such a way that { 7 q} is a non-negative P-supermartingale. By 


the equations (2.2) and (2.3), it is straightforward to see that 


PtT = 


Rqt + b\{T)A^ 

Pot + b^A^ ' 


(i) 


n = - 


Rot + b\ (t)A 
P ot + b 1 (t)A[ iy 


(2.15) 


where the “dot-notation” means differentiation with respect to time t. 
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Let us return to the modelling of rational multi-curve term structures and in particular 


to the definition of the (forward) LIBOR process. Putting equations (2.6) and (2.1) in 


relation, we see that the model (2.13) naturally offers itself as a model for the LIBOR 


process (2.6) in the considered setup. Since (2.13) satisfies (2.1) by construction, so does 


the LIBOR model 




L(0‘, Tj_i,Ti) + b 2 (T l _ 1 ,T i )A ( i 2) + b 3 (^-i,r ; ;)4 3) 

Pot + b^A^ 


(2.16) 


satisfy the martingale equation (2.6) and in particular (2.10) for t < Tj_ i. In Macrina (2014|) 
a method based on the use of weighted heat kernels is provided for the explicit construction 
of the M-martingales {^4[^}i=i,2 and thus in turn for explicit LIBOR processes. The method 
allows for the development of LIBOR processes, which, if circumstances in financial markets 
require it, by construction take positive values at all times. 


2.3 Bottom-up risk-neutral approach 

Since we also deal with counterparty-risk valuation adjustments, we present another scheme 
for the construction of the LIBOR models, which we call “bottom-up risk-neutral ap¬ 
proach”. As the name suggest, we model the multi-curve term structure by making use 
of the risk-neutral measure (via the auxiliary measure M) while the connection to the 
P-dynamics of prices can be reintroduced at a later stage, which is important for the cal¬ 
culation of risk exposures and their management. “Bottom-up” refers to the fact that the 
short interest rate will be modelled first, then followed by the discount bond price and LI¬ 


BOR processes. Similarly, in Section 6.1, the hazard rate processes for contractual default 
will be modelled first, and thereafter the price processes of counterparty risky assets will be 
derived thereof. We utilise the notation E[... \Pt\ = E*[...]. In the bottom-up setting, we 


directly model the short risk-free rate {ry} in the manner of the right-hand side in (2.15), 


i.e. 


n = 


c\{t) + bijt)^ 
ci(t) + bi(t)A 


d) : 


(2.17) 


by postulating (i) non-increasing deterministic functions b\(t) and c\(t) with ci(0) = 1 
(later c\(t) will be seen to coincide with Pot); and (ii) an ({P)},M)-martingale {A^} with 


Aq 1 ^ = 0 such that 


is a positive ((p}, 


ht = ci(t) + bi(t)A\ 
)-supermartingale for all t > 0. 


(i) 


(2.18) 


Example 2.1. Let A) ±; = S)' 1 — 1, where {SjW} is a positive M-martingale with Sg 1 * = 1; 


for example exponential Levy martingales. The supermartingale (2.18) is positive for any 
given t if 0 < b\(t) < c\{t). 


Associated with the supermartingale (2.18), we characterise the (risk-neutral) pricing mea¬ 
sure Q by the M-density process {fit}o<t<T, given by 


/A = 




= £ 


r fciWdA^ 

'o c\(t) + &i(t)A^ ; 


(2.19) 
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which is taken to be a positive ({J-)}, M)-martingale. Furthermore, we denote by Dt = 
exp Jq r s ds'j the discount factor associated with the risk-neutral measure Q. 

Lemma 2.1. h = Dt Ht- 

Proof. The Ito semimartingale formula applied to ip(t, A = ln(ci (t) + bi(t)A[^) = In (lit) 
and to In (DtHt) gives the following relations: 


din ^ci(t) + bi^A^ 1 ^ = — rtdt + 


b\{t)dA 


(i) 


bl{t)d[A^\AW] 


ci(t) + biitfA^ 1 } 2(ci(t) -I- bitPjA^}) 2 

6i(t)A4 1} 


+ d^ Ain (ci(f) + bi(t)Af’) - 


S<t 


ci(t) + bi{t)A^} / 


( 2 . 20 ) 


where (2.17) was used in the first line, and 


dln(DtfJ,t) = dlnDt + dlnfit 




A ^ 


= -r t dt + 


2 (Ht-) 2 
bi(t)dAf 1 ' 1 

ci{t) + bi(t)A^} 2(ci(t) + 6i(t)y4^) 2 


Ht- 

b\{t)d[AA),AA)f t 


( Aln(^ t 


S<t 


where 


A In (fit) = In ( -p- ) = In ( 1 + 


d(t) + b 1 {t)A^ 1 } / 
= Ain (d(t) + bityAj 1 ^ . 


= In 


M*)AA[ 1} 

ci(t) + bi(t)A^ 

d(t) + b 1 (t)A^ 1) \ 

^d(t) + b 1 (t)Aj 1) J 


( 2 . 21 ) 


Therefore, dln(ht) = dln(Th^). Moreover, h o = Dqho = 1. Hence h t = D t Ht- ■ 

It then follows that the price process of the OIS discount bond with maturity T can be 
expressed by 


PtT = 


D r -p 
~D~t 


1 


D t Ht 


E M [D t Ht | Ft] = Ef 


hx 


dm+hWA™ 

ci(t) + bi(t)A 


(i) 


( 2 . 22 ) 


for 0 < f < T, Thus, the process {ht} plays the role of the pricing kernel associated with 
the OIS market under the measure M. In particular, we note that d(t) = Pot for t £ [0,T] 
and r’t = — (dr In PtT)\x=t — 0- A construction inspired by the above formula for the OIS 
bond leads to the rational model for the LIBOR prevailing over the interval [Tj_i, Tf). The 
Ft,_ x -measurable spot LIBOR rate L(Ti\ Tj_i, Tf) is modelled in terms of } and, in 
this paper, at most two other M-martingales {A) '} and {A t ; } evaluated at Tj_j: 


L(0-,Ti_i,Ti) + b 2 {T i _ 1 ,T i )A^_ i + b 3 {T i _ 1 ,T i )A^>_ 
PoTi + h{Tf)A^ 


,( 3 ) 


T(T); Ti— i, Tf) — 


(2.23) 
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The (forward) LIBOR process is then defined by an application of the risk-neutral valuation 
formula (which is equivalent to the pricing formula (2.1) under P) as follows. For t < Tj_i 
we let 


LfaTi-uTi) = ±-E®[D Ti L(T i -T i _ 1 ,T i )}=Ef 

Dt 


D Tj VTj 

A IH 


L{Ti\Ti_i,Ti) 


= E5 


ht 


(2.24) 

(2.25) 


and thus, by applying (2.18) and (2.23), 


L(t;Tj_i,Tj) — 


L{ 0; Ti i, Ti) + 


Pot + bi(t)A 


(i) 


(2.26) 


Hence, we recover the same model (and expression) as in (2.16). The LIBOR models (2.26) 
(or (2. 16|)) are compatible with an HJM multi-curve setup where, in the spirit of Heath, 


Jarrow and Morton (1992), the initial term structures Pot* and L(0;T,_i,Tj) are fitted by 
construction. 

Example 2.2. Let A^ = — 1, where sj^ is a positive M-martingale with S^' 1 = 1. 

For example, one could consider a unit-initialised exponential Levy martingale defined in 
terms of a function of an M-Levy process for i = 2,3. Such a construction produces 

non-negative LIBOR rates if 


0 < f>2(A—i, Ti) + bz(Ti_i,Ti) < L( 0; Ti-i, T*). 


(2.27) 


If this condition is not satisfied, then the LIBOR model may be viewed as a shifted model, 
in which the LIBOR rates may become negative with positive probability. For different 


kinds of shifts used in the multi-curve term structure literature we refer to, e.g., Mercurio 
|(2010a ) or Moreni and Pallavicini (2014). 


3 Clean valuation 


The next questions we address are centred around the pricing of LIBOR derivatives and 
their calibration to market data, especially LIBOR swaptions, which are the most liquidly 
traded (nonlinear) interest rate derivatives. Since market data typically reflect prices of 
fully collaterallised transactions, which are funded at a remuneration rate of the collateral 
that is best proxied by the OIS rate, we consider in this section, in the perspective of model 
calibration, clean valuation ignoring counterparty risk and assume funding at the rate rt- 


An interest rate swap (see, e.g., Brigo and Mercurio (2006)) is an agreement between 
two counterparties, where one stream of future interest payments is exchanged for another 
based on a specified nominal amount N. A popular interest rate swap is the exchange of 
a fixed rate (contractual swap spread) against the LIBOR at the end of successive time 
intervals [T)_i, T,;j of length 5. Such a swap can also be viewed as a collection of n forward 
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rate agreements. The swap price Swt at time t < Tq is given by the following model- 
independent formula: 

n 

Sw t = Nd^mt-T^Ti) - KP m }. 

1=1 

A swaption is an option between two parties to enter a swap at the expiry date T\ (the 
maturity date of the option). Its price at time t < T*. is given by the following M-pricing 
formula: 


Swn tTk = -r^E M [h Tk (Sw Tk ) + \J : 't\ 
Ht 


h t 


h Tk E [HT^T^Tt) - KP TkTi ] 




Ft 


Nb 


Pot + A t 


(i) 


E k 


m 

( E [^(0;T i _i,r i ) + 6 2 (T i _i,T i )Ag ) + 6 3 (iT.-i,T i )A. 


(3) 

T k 


i=k-\-1 


-KiPw+hWAg 


Ft 


(3.28) 


using the formulae (2.22) and (2.26) for Pr k T- and L(Tj.; Tj_i, Tj). In particular, the swap- 

- - • 1 (£) (i) 

tion prices at time t = 0 can be rewritten by use of A\’ = S t — 1 so that 


Swn 0Tk 


AW 

N5F} 


(c 2 A, 




C2S%+C 3 ^ 


Cij4 A + Co) 

+ c 0 ) , 


(3.29) 


where 

mm m 

E & 2 (Tj_i,Tj), c 3 = E & 3(7i-i,T)), C] = I< E 

i=fc+l i=fc+l z=fc+l 

m 

E [-k(0> Ti-l>Ti) - KPoTi ], Co = Co + Cl — c 2 — C3. 

2=fc+l 

As we will see in several instance of interest, these expectations can be computed efficiently 
with high accuracy by various numerical schemes. 

Remark 3.2. The advantages of modelling the LIBOR process {L(f; T,- 1 . Tj )} by a rational 
function of which denominator is the discount factor (pricing kernel) associated with the 
employed pricing measure (in this case M) are: (i) The rational form of {L(t; T)_i, Tj)} and 
also of {PtTi} produces, when multiplied with the discount factor {ht}, a linear expression in 
the M-martingale drivers {A^}. This is in contrast to other akin pricing formulae in which 
the factors appear as sums of exponentials, see e.g. |Crepey, Grbac, Ngor and Skovmand 
(2014]), Equation (33). (ii) The dependence structure between the LIBOR process and 
the OIS discount factor {ht }—or the pricing kernel {-/r t} under the P-measure—is clear- 
cut. The numerator of {L(f; Tj_i, Tj)} is driven only by idiosyncratic stochastic factors 
that influence the dynamics of the LIBOR process. We may call such drivers the “LIBOR 


C2 = 

co = 
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risk factors”. Dependence on the “OIS risk factors”, in our model example {Aj 1 }, is 
produced solely by the denominator of the LIBOR process, (iii) Usually, the FRA process 
Kt = L(t; Ti -\, T,')/PtT. t is modelled directly and more commonly applied to develop multi¬ 
curve frameworks. With such models, however, it is not guaranteed that simple pricing 


formulae like (3.28) can be derived. We think that the “codebook” (2.6), and (2.26) in 


the considered example, is more suitable for the development of consistent, flexible and 
tractable multi-curve models. 


3.1 Univariate Fourier pricing 

Since in current markets there are no liquidly-traded OIS derivatives and hence no useful 
data is available, a pragmatic simplification is to assume deterministic OIS rates rt■ That 
is to say = 0 , and hence b\(t) plays no role either, so that it can be assumed equal to 

/ o\ _ 

zero. Furthermore, for a start, we assume A\ = 0 and 63 (t) = 0, and (3.29) simplifies to 


SwriQT k = N5 E 


M 


(c 2 A 


( 2 ) 


+ Co 


= N5W 


C 2 Sg } +C 0 ) 


where here Co = co — c 2 . For Co > 0 the price is simply SwriQT k = N5cq. For Co < 0, and in 
the case of an exponential-Levy martingale model with 

s (2) = e xV-tM L 

(o) 

where {X) '} is a Levy process with cumulant i/j 2 such that 


E 




( 2 ) 


= exp [tij) 2 {z)\ , 


we have 


where 


N5 

Swn QTk = — 


1—iv—R 


0 


M^\r + iv) 


2ir J K (R + iv)(R + iv — 1) 

M^(z) = e Tfe ^ 2 (A+ j: ( lri (c2)-i/’2(i)) 
T k ^ ' 


dv. 


(3.30) 

(3.31) 


(2) 

and R is an arbitrary constant ensuring finiteness of Mrp k (R + iv) for v E M. For details 


concerning (3.31), we refer to, e.g., Eberlein, Glau and Papapantoleon (2010). 


3.2 One-factor lognormal model 

In the event that {A^} = {Aj 3 )} = 0 and {Ais of the form 


A{ 2) = exp ( a 2 x[ 2) - ) - 1 , 


(3.32) 


(2) 

where {A) '} is a standard Brownian motion and a 2 is a real constant, it follows from 
simple calculations that the swaption price is given, for cq = cq — c 2 , by 


Swn m =N5W‘ 


{c 2 A^)+c^j 


atz( ^ (\a 2 2 T -\n(c 0 /c 2 )\ 

=m ^ y - -j= -j + c 0 d> 


-\cl\T - ln(c 0 /c 2 )\ 


a 2 VT 


(3.33) 

(3.34) 


where <h(x) is the standard normal distribution function. 
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3.3 Two-factor lognormal model 


We return to the price formula (3.29) and consider the case where the martingales {nj*' 1 } 
are given, for i = 1,2,3, by 


A\ l) = exp 




(3.35) 


for real constants a* and standard Brownian motions {X^} = and {X^} with 

correlation p. Then it follows that 


Swn 0Tk = E M 


^ C2 g' Y VTfca2 — \a%T k _|_ 'JTk a 3 — ^af i T k _ ^^Yy/T^ai 


2 a l 1 k 


+ Cq 


(3.36) 


where X ~ jV(0,1), Y ~ jV(0,1), (X|H) = y ~ M(py , (1 — p 2 )). Hence, 

/ oo /*oo 

/ (c 2 e a: ' /7 ^“ 2_ 2 a 2 Tfc - J\(y)) + /(x|?/)/(y)dxdy 

-OO J — OO 

= [ ([ (c 2 e x%/Tka2 ~^ Tk - K(y)) + f(x\y)dx\ f{y)dy 

■JK(y)>0 Vi —oo / 

+ [ ([ (c 2 e Xy/Tka2 ~^ Tk - K(y)) + f(x\y)dx] f(y)dy, 

JK(y )<0 \J—oo / 


where 


if(y) = Cl (e ai ^-^ a ? Tfc - 1) - C3 ( e ° 3 ^-|a§T fc _ j) _ C0; 


/(y) = 
/(*|y) 


V2tt 


\/2vr(l - p 2 ) 


-(x-py ) 2 

> 2(1-P 2 ) 


This expression can be simplified further to obtain 


SwriQT k 


’ K{y)> o 


C2e 


a2 ^ ra+ i a 2 Tfc (i-p 2 )$ f PV + Q 2 V^fe(l ~ P 2 ) + ln(c 2 ) - \alT k - K(y) 


-K(y )$ 


py + ln(c 2 ) - ga 2 Tfc - if (y) 


V 1 - p 2 


f(y)dy 


+ 


lK(y)<0 


c 2 e 


a.2VTkP(y-bo'2VTkP _ 


if(y)) f(y)dy- 


The calculation of the swaption price is then reduced to calculating two one-dimensional 
integrals. Since the regions of integration are not explicitly known, one has to numerically 
solve for the roots of K(y), which may have up to two roots. Nevertheless a full swaption 
smile can be calculated in a small fraction of a second by means of this formula. 
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4 Calibration 


The counterparty-risk valuation adjustments, abbreviated by XVAs (CVA, DVA, LVA, etc.), 
can be viewed as long-term options on the underlying contracts. For their computation, 
the effects by the volatility smile and term structure matter. Furthermore, for the planned 
XVA computations of the multi-curve products in Section [ 6 j it is necessary to calibrate the 
proposed pricing model to financial instruments with underlying tenors of 5 = 3m and 5 = 
6 m (the most liquid tenors). Similar to Crepey, Grbac, Ngor and Skovmand (2014), we make 
use of the following EUR market Bloomberg data of January 4, 2011 to calibrate our model: 
EONIA, three-month EURIBOR and six-month EURIBOR initial term structures on the 
one hand, and three-month and six-month tenor swaptions on the other. As in the HJM 
framework of Crepey, Grbac, Ngor and Skovmand (2014), to which the reader is referred for 
more detail in this regard, the initial term structures are fitted by construction in our setup. 
Regarding swaption calibration, at first, we calibrate the non-maturity/tenor-dependent 
parameters to the swaption smile for the 9x1 years swaption with a three-month tenor 
underlying. The market smile corresponds to a vector of strikes [—200,—100,—50,—25, 
0, 25, 50,100, 200] bps around the underlying swap spread. Then, we make use of at-the- 
money swaptions on three and six-month tenor swaps all terminating at exactly ten years, 
but with maturities from one to nine years. This co-terminal procedure is chosen with a 
view towards the XVA application in Section [ 6 j where a basis swap with a ten-year terminal 
date is considered. 

In particular, in a single factor {A\ ; } setting: 


( 2 ) 

1 . First, we calibrate the parameters of the driving martingale {A\ to the smile of the 
9x1 years swaption with tenor 5 = 3m. This part of the calibration procedure gives 
us also the values of 62 ( 9 , 9 . 25 ), 62 ( 9 . 25 , 9 . 5 ), 62 ( 9 . 5 , 9 . 75 ) and 62 ( 9 . 75 ,10), which we 
assume to be equal. 


2. Next, we consider the co-terminal, A x (10 — A), ATM swaptions with A = 1,2,..., 
9 years. These are available written on the three and six-month rates. We calibrate 
the remaining values of 62 one maturity at a time, going backwards and starting with 
the 8x2 years for the three-month tenor and with the 9x1 years for the six-month 
tenor. This is done assuming that the parameters are piecewise constant such that 
62 (T, T + 0.25) = 6 2 (T + 0.25, T + 0.5) = 6 2 (T + 0.5, T + 0.75) = 6 2 (T + 0.75, T + 1) 
for each T = 0,1,..., 8 and that 62 ( 1 ", T + 0.5) = bs(T + 0.5, T + 1) hold for each 
T = 0,l,...,9. 


4.1 Calibration of the one-factor lognormal model 


In the one-factor lognormal specification of Section |3.2[ we calibrate the parameter a 2 
and 6 = 62 ( 9 , 9.25) = 62 ( 9 . 25 , 9 . 5 ) = 62 ( 9 . 5 , 9 . 75 ) = 62 ( 9 . 75 , 10 ) with Matlab utilising the 
procedure “lsqnonlin” based on the pricing formula (3.33) (if cq < 0, otherwise SwnoT k = 
N5cq). This calibration yields: 


a 2 = 0.0537,6 = 0.1107. 
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Forcing positivity of the underlying LIBOR rates means, in this particular case, restricting 
6 < L(0;9.75,10) = 0.0328 (cf. (2.27)). The constrained calibration yields: 


a 2 = 0.1864,6 = 0.0328. 


The two resulting smiles can be found in Figure [lj where we can see that the unconstrained 
model achieves a reasonably good calibration. However, enforcing positivity is highly re¬ 
strictive since the Gaussian model, in this setting, cannot produce a downward sloping 
smile. 



Strike 

Figure 1: Lognormal one-factor calibration 


Next we calibrate the 6 2 parameters to the ATM swaption term structures of 3 months 
and 6 months tenors. The results are shown in Figure [2j When positivity is not enforced 
the model can be calibrated with no error to the market quotes of the ATM co-terminal 
swaptions. However, one can see from the figure that the positivity constraint does not 
allow the 6 2 function to take the necessary values, and thus a very poor fit to the data is 
obtained, in particular for shorter maturities. 

With this in mind the natural question is whether the positivity constraint is too re¬ 
strictive. Informal discussions with market participants reveal that positive probability for 
negative rates is not such a critical issue for a model. As long as the probability mass for 
negative values is not substantial, it is a feature that can be lived with. Indeed assign¬ 
ing a small probability to this event may even be realistic 0 In order to investigate the 
significance of the negative rates and spreads mentioned in remark 2.1, we calculate lower 
quantiles for spot rates as well as the spot spread for the model calibrated without the 
positivity constraint. As Figure [3] shows, the lower quantiles for the rates are of no concern. 
Indeed it can hardly be considered pathological that rates will be below -14 basis points 
with 1% probability on a three year time horizon. Similarly, with regard to the spot spread, 
the lower quantile is in fact positive for all time horizons. Further calculations reveal that 
the probabilities of the eight year spot spread being negative is 1.1 x 10~' 5 and the nine 
year is 0.008 - which again can hardly be deemed pathologically high. 


J As with EONIA since the end of 2014 or Swiss rates in the crisis. 
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We find that the model performs surprisingly well despite the parsimony of a one-factor 
lognormal setup. While positivity of rates and spreads are not achieved, the model assigns 
only small probabilities to the negatives. However, the ability of fitting the smile with such 
a parsimonious model is not satisfactory (cf. Fig. [I]), which is our motivation for the next 
specification. 





Expiry in years T 

Figure 2: One-Factor Lognormal calibration. (Left) Fit to ATM swaption implied volatility 
term structures. (Right) Calibrated values of the 62 parameters. (Top) 5 = 3m. (Bottom) 
5 = 6 m. 
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Figure 3: One-Factor Lognormal calibration. 1% lower quantiles 
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4.2 Calibration of exponential normal inverse Gaussian model 

( 2 ) 

The one-factor model, which is driven by a Gaussian factor {A) is able to capture the 
level of the volatility smile. Nevertheless, the model implied skew is slightly different from 
the market skew. To overcome this issue, we now consider a one-factor model driven by a 
richer family of Levy processes. The process {A\ '} is now assumed to be the exponential 
normal inverse Gaussian (NIG) M-martingale 

A? = exp (x f (2) - ^(1)) - 1, (4-37) 


( 2 ) - 

where {X) '} is an M-NIG-process with cumulant ip(z), see (3.30), expressed in terms of 

the parametrisatiorj^] (u, 9 , a) from Cont and Tankov (2003) as 

tj){z ) = —v (\/v 2 — 2 z6 — z 2 a ' 2 — iA , (4.38) 


where u, a > 0 and 8 E M. The parameters that need to be calibrated at first are v, 9 , a 
and b = 62(9,9.25) = 6o(9.25,9.5) = 62 ( 9 .5, 9.75) = 62 ( 9 . 75 ,10). After the calibration, we 
obtain 


6 = 0.0431, v = 0.2498, 8 = -0.0242, <r = 0.1584. 


Imposing 6 < L{ 0; 9.75,10) = 0.0328 to get positive rates we obtain instead 


6 = 0.0291, 1 / = 0.1354, 8 = -0.0802, cr = 0.3048. 

The two fits are plotted in Figure[4j Here, imposing positivity comes at a much smaller cost 
when compared to the one factor Gaussian case. The NIG process has a richer structure 
(more parametric freedom) and therefore is able to compensate for an imposed smaller level 
of the parameter 62 . 

2 The 



parametrisation is recovered by setting (l — 0. 


1 

a = — 
a 


a? +Vi 


8 = 


15 















Strike 

Figure 4: Exponential-NIG calibration 


We continue with the second part of the calibration of which results are found in Figure 
[5] Here we see that enforcing positivity may have a small effect on the smile but it means 
that the volatility structure cannot be made to match swaptions with maturity smaller 
than 7 years. Thus, enforcing positivity in this model produces limitations which we wish 
to avoid. In Figure [6] we plot lower quantiles for the rates and spreads as for the one-factor 
lognormal model. While spot spreads remain positive, the levels do not, and, as shown, the 
model assigns an unrealistically high probability mass to negative values. In fact the model 
assigns a 1% probability to rates falling below -12% within 2 years! Thus, the one-factor 
exponential-NIG model loses much of its appeal for it cannot, in a realistic manner, be 
made to fit long-term smiles and shorter-term ATM volatilities. 


4.3 Calibration of a two-factor lognormal model 


The necessity to produce a better fit to the smile than what can be achieved with the one- 
factor Gaussian model, while maintaining positive rates and spreads, leads us to proposing 
the two-factor specification presented in Section 3.3 This model is heavily parametrised 


and the parameters at hand are not all identified by the considered data. We therefore fix 
the following parameters: 


ai = 1, a3 = 1.6, (4.39) 

6 3 (T ,T + 0.25) = 0.15L(0;T;T + 0.25), Te [9,9.75], (4.40) 

b 2 (T, T + 0.25) = 0.55L(0,T;T + 0.25), Te [0,8.75]. (4.41) 

We assume that b\ is constant, i.e. b\ = b\{T ) for T e [0,10], and that 63 , outside 
of the region defined above, is piecewise constant such that 63 (T, T + 0.25) = 63 (T + 

0.25, T + 0.5) = b 3 {T + 0.5, T + 0.75) = b 3 (T + 0.75, T + 1) for each T = 0,1..., 8 and 


b 3 (T, T + 0.5) = b 3 (T + 0.5, T + 1) holds for each T = 0,1..., 9. We furthermore assume 
that b 2 (T,T + 0.5) = b 2 (T,T + 0.25), T € [0,9.5]. These somewhat ad hoc choices are 
done with a view towards b 2 and b 3 being fairly smooth functions of time. We herewith 
apply a slightly altered procedure to calibrate the remaining parameters if compared to the 
scheme utilised for the one-factor models. 
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Figure 5: Exponential-NIG calibration. (Left) Fit to ATM swaption implied volatility term 
structures. (Right) Calibrated values of the 62 parameters. (Top) 5 = 3m. (Bottom) 
5 = 6m. 
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Figure 6 : Exponential-NIG calibration calibration. 1% lower quantiles. 
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1. We first calibrate to the smile of the 9x1 years swaption which gives us the param¬ 
eters 0,2 , Pi the assumed constant value of b\, and 62 ( 9 , 9 - 25 ) to 62 ( 9 .75,10) which are 
assumed equal to a constant 6. Similar to the exponential-NIG model, we make use 
of four parameters in total to fit the smile. 

2. The remaining 62 parameters are determined a priori , so what remains is to cali¬ 
brate the values of 63 . The three-month tenor values bs(T,T + 0.25) for T e [0,8.75] 
are calibrated to ATM, co-terminal swaptions starting from the 8x2 years and then 
continuing backwards to the 1x9 years instruments. For the six-month tenor prod¬ 
ucts, we calibrate bs(T,T + 0.5) for T G [0,9.5] starting with 9x1 years and proceed 
backwards. 

These are the values we obtain from the first calibration phase: 61 = 0.2434, 6 = 

0.02, 02 = 0.1888, p = 0.9530. The corresponding fit is plotted in the upper left quad¬ 

rant of Figure [7J In order to check the robustness of the calibrated fit through time, we 




Strike Strike 




Strike Strike 


Figure 7: Lognormal two-factor calibration. 


also calibrate to three alternative dates. The quality of the fit appears quite satisfactory 
and comparable to the exponential-NIG model. For all four dates the calibration is done 
enforcing the positivity condition b 2 {T,T + 0.25) + bs(T,T + 0.25) < L(0;T, T + 0.25). 
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However, the procedure yields the exact same parameters even if the constraint is relaxed. 
We thus conclude that a better calibration appears not to be possible for these datasets 
by allowing negative rates. Note that it is only for our first data set that the calibrated 
correlation p is as high as 0.9530. In the other three cases we have p = 0.4118, p = 0.3964, 
and p = 0.2461. Figure [8] shows the parameters 62 and 63 obtained at the second phase 
of the calibration to the data of 4 January 2011. As with the previous model (cf. the left 
graphs of Figures [ 2 ] and [5]), the volatilities are matched to market data without any error. 




Figure 8 : Two-factor lognormal calibration. (Left) Parameter values fitted to three-month 
ATM swaption implied volatility term structures. (Right) Parameter values fitted to six- 
month ATM swaption implied volatility term structures. 


We add here that, although not visible from the graphs, the calibrated parameters 


satisfy the LIBOR spread positivity discussed in Remark 2.1 


In conclusion, we find that the two-factor log-normal has the ability to fit the swaption 
smile very well, it can be controlled to generate positive rates and positive spreads, and it is 
tractable with numerically-efficient closed-form expressions for the swaption prices. Given 
these desirable properties, we discard the one-factor models and retain the two-factor log¬ 
normal model for all the analyses in the remaining part of the paper. 


5 Basis swap 

In this section, we prepare the ground for counterparty-risk analysis, which we shall treat 
in detail in Section [ 6 j A typical multi-curve financial product, i.e. one that significantly 
manifests the difference between single-curve and a multi-curve discounting, is the so-called 
basis swap. Such an instrument uconsists of exchanging two streams of floating payments 
based on a nominal cash amount N or, more generally, a floating leg against another 
floating leg plus a fixed leg. In the classical single-curve setup, the value of a basis swap 
(without fixed leg) is zero throughout its life. Since the onset of the financial crisis in 2007, 
markets quote positive basis swap spreads that have to be added to the smaller tenor leg, 
which is clear evidence that LIBOR is no longer accepted as an interest rate free of credit 
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or liquidity risk. We consider a basis swap with a duration of ten years where payments 
based on LIBOR of six-month tenor are exchanged against payments based on LIBOR of 
three-month tenor plus a fixed spread. The two payment streams start and end at the same 
times To = Tg 1 = Tg, T = T’ = T 2 9 . The value at time t of the basis swap with spread K 
is given by 


n i 


ri2 


BS t = N PT <5f m T(t; Tl 1: T}) - E Tj^, T 2 ) + KP tr? 

\*=1 3 =1 

for t < To. After the swap has begun, i.e. for To < t < T, the value is given by 

( ni 

5f™L(Tl_ 1 ;Tl_ 1 ,Tl)+ E i’ T i) 

i=it+1 

n 2 

~(j- 1 (T/ t _ i; Tj t _i , T ? t ) + KP tT 2 ^) — E S j m ( L (t; T 2 _1, T 2 ) + KP t 

j=jt+1 

where T- f (respectively T 2 ) denotes the smallest Tl (respectively Tf) that is strictly greater 
than t. The spread K is chosen to be the fair basis swap spread at To so that the basis 
swap has value zero at inception. We have 



K = 


Efi, Sf'm-, TU,T}) - Eji, sf “L(T„; 

i ^-Pt.t? 


The price processes on which the numerical illustration in Figure [9] have been obtained 
was simulated by applying the calibrated two-factor lognormal model developed in Section 
4.3 The basis swap is assumed to have a notional cash amount N = 100 and maturity 


T = 10 years. In the two-factor lognormal setup, the basis swap spread at time t = 0 
is K = 12 basis points, which is added to the three-month leg so that the basis swap is 
incepted at par. The t = 0 value of both legs is then equal to EUR 27.96. The resulting risk 
exposure, in the sense of the expectation and quantiles of the corresponding price process at 
each point in time, is shown in the left graphs of Figure [9j where the right plots correspond 
to the P exposure discussed in Section |5Tj Due to the discrete coupon payments, there are 
two distinct patterns of the price process exposure, most clearly visible at times preceding 
payments of the six-month tenor coupons for the first one and at times preceding payments 
of the three-month tenor coupons without the payments of the six-month tenor coupons 
for the second one. We show the exposures at such respective dates on the upper and lower 
plots in Figure [9| 
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Basis swap positive exposure under M-measure Basis swap positive exposure under P-measure 



Figure 9: Exposures of a basis swap (price process with mean and quantiles) in the cali¬ 
brated two-factor Gaussian model. (Top) Exposure of the basis swap at t = 5m, 11m, etc. 
(Bottom) Exposure of the basis swap price at t = 2m, 8m, 14 m, etc. (Left) Exposure under 
the M-measure. (Right) Exposure under the P-measure with the prediction that LIBOR 
rate L(10.75y; 10.75y, lly) will be either 2% with probability p = 0.7 or 5% with proba¬ 
bility 1 — p = 0.3. 

5.1 Levy random bridges 

The basis swap exposures in Figure [9] are computed under the auxiliary M-measure. The 
XVAs that are computed in later sections are derived from these M-exposures. However, 
exposures are also needed for risk management and as such need to be evaluated under the 
real-world measure P. This means that a measure change from M to P needs to be defined, 
which requires some thoughts as to what features of a price dynamics under P one might like 
to capture through a specific type of measure change and hence by the induced P-model. In 
other words, we design a measure change so as to induce a particular stochastic behaviour 
of the {At} processes under P, and in particular of the underlying Markov processes {Xt} 
driving them. 

A special case we consider in what follows is where {A 7 *} is a Levy process under M, 
while it adopts the law of a corresponding (possibly multivariate, componentwise) Levy 
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random bridge (LRB) under P. Several explicit asset price models driven by LRBs have 


been developed in Macrina (2014). The LRB-driven rational pricing models have a finite 


time horizon. The LRB is characterised, apart from the type of underlying Levy process, 
by the terminal P-marginal distribution to which it is pinned at a fixed time horizon U. 
The terminal distribution can be arbitrarily chosen, but its specification influences the 
behaviour of the LRB as time approaches U. In turn, the properties of a specified LRB 
influence the behaviour of {At} and hence the dynamics of the considered price process. 
We see an advantage in having the freedom of specifying the P-distribution of the factor 
process at some fixed future date. This way, we can implement experts’ opinions (e.g. 
personal beliefs based on some expert analysis) in the P-dynamics of the price process as 
to what level, say, an interest rate (e.g. OIS, LIBOR) is likely to be centred around at a 
fixed future date. 


The recipe for the construction of an LRB can be found in Hoyle, Hughston and Macrina 
(2011), Definition 3.1, which is extended for the development of a multivariate LRB in 
|Macrina (2014). LRBs have the property, as shown in Proposition 3.7 of Hoyle, Hughston 


and Macrina (2011), that there exists a measure change to an auxiliary measure with respect 


to which the LRB has the law of the constituting Levy process. That is, we suppose the 
auxiliary measure is M and we have an LRB {W}o<t<l/ defined on the finite time interval 
[0,17] where U is fixed. Under M and on [0,17), {7Q} has the law of the underlying 
Levy process. To illustrate further, let us assume a univariate LRB; the analogous measure 


change for multivariate LRBs is given in Macrina (2014). Under P, which stands in relation 
with M via the measure change 


Vt = 


dP 

dM 


Tt 


I 

7m 


fu-t(z - X t ) 
fu(z) 


v(&z), t < U, 


(5.42) 


where ft(x) is the density function of the underlying Levy process for all t E (0,17] and v is 
the P-marginal law of the LRB at the terminal date U, the process {Xt} is an LRB (note 
that the change of measure is singular at U ). 

Now, returning to the calibrated two-factor lognormal model of Section |4.3[ but similarly 
also to the other models in Section Ul we may model the drivers {xf 1 ^} = {xj 3 '*} and 
{A') by two dependent Brownian random bridges under P. The computed M exposures 
in Figure [9] thus need to be re-weighted by the corresponding amount ijt in order to obtain 
the P-exposures of the basis swap. Since here we employ LRBs, we have the opportunity 
to include an expert opinion through the LRB marginals v as to what level one believes 
the interest rates will tend to by time U. The re-weighted P-exposures of the basis swap 
are plotted in the graphs of the right-hand side of Figure [9] The maximum of the upper 
quantile curves shown in the graphs is known as the potential future exposure (PFE) at 
the level 97.59c[f} 

Hence, we now have the means to propose a risk-neutral model that can be calibrated to 
option data, and which after an explicit measure change can be applied for risk management 
purposes while offering a way to incorporate economic views in the dynamic of asset prices. 

3 In practice, people rather consider the expected positive exposure (expectation of the positive part of 
the price rather than the price) in the PFE computation, but the methodology is the same. 
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Recalling (2.19) and (5.42), the Q-to-P measure change is obtained by 


dP 


Ft vt 


(5.43) 


and the pricing formula for financial assets ( 2 . 1 ) may be utilised under the various measures 
as follows: 


StT = [Dt Stt | F t ] = [D t v t Stt | Ft] = [h T S TT \ F t ] 

D t D t v t h t 


Vt -E p 


D t v t 


Dt vt 
Vt 


Stt 


Ft 


= —E f [7tt Stt \ Ft], 

Tt 


(5.44) 


for 0 < t < T < U (since we consider price models driven by LRBs). It follows that the 
pricing kernel is given by 77 = DtVtr]^ 1 = rj^ 1 ht- Measure changes from a risk-neutral 
to the real-world probability measure are discussed for similar applications also elsewhere. 
For a recent study in this area of research, we refer to, e.g., 


Hull, Sokol and White (2014 


6 Adjustments 

So far we have focused on so-called “clean computations”, i.e. ignoring counterparty risk 
and assuming that funding is obtained at the risk-free OIS rate. In reality, contractually 
specified counterparties at the ends of a financial agreement may default, and funding to 
enter or honour a financial agreement may come at a higher cost than at OIS rate. Thus, 
various valuation adjustments need to be included in the pricing of a financial position. 
The price of a counterparty-risky financial contract is computed as the difference between 
the clean price, as in Section [3j and an adjustment accounting for counterparty risk and 
funding costs. 


6.1 Rational credit model 


As we shall see below, in addition to their use for the computation of PFE, the exposures in 
Section [5] can be used to compute various adjustments: CVA (credit valuation adjustment), 
DYA (debt valuation adjustment) and LVA (liquidity-funding valuation adjustment). With 


this goal in mind, we equip the bottom-up construction in Section 2.3, the notation of which 


is used henceforth, with a credit component in the following manner. 

-Mr*—i, 2 ,...,n, w j 1 j c j 1 are assumed to be ({Ft}, M)-Markov processes. For 


We consider {X [ t ; } 0<f 
any multi-index (*i,we write F t 


(*!,■••,*d) 


= Vi=i,...,d^ 


■x(h) 


The (market) filtration 


{Ft} is given by {jr(b—> n )}. p or the application in the present section, we fix n = 6 . 
The Markov processes {A'^} = {A'( 3 )} and {X^} are utilised to drive the OIS and 
LIBOR models as described in Section |2.3l in particular the zero-initialised ({.Tt},M)- 


martingales 1 > 2 ’ 3 . The Markov processes i = 4,5,6, which are assumed to 

be M-independent betweeen them and of the Markov processes i = 1,2,3, are applied to 


model {.Til-adapted processes { 7 ^}® 4,5,6 defined by 


M 

Tt = 


Ci(t ) + bi(t)A 


(*) 


Ci(t) + bi(t)A 


M 


(6.45) 
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where bi{t) and Cj(f), with Cj(0) = 1, are non-increasing deterministic functions, and where 
{i 4^}* =4 ’ 5 ’ 6 are zero-initialised ({7*}, M)-martingales of the form A(t,X^). Comparing 
with ( [2l7| ), we see that (|6.45[) is modelled in the same way as the OIS rate (j2.17[), non¬ 


negative in particular, as an intensity should be (see Remark 6.3). 


In line with the “bottom-up” construction in Section (2.3), we now introduce a density 


({71}, M)-martingale {^t^t}o<t<T that induces a measure change from M to the risk-neutral 


measure 


Tt 


= TtPt (0 < t < T) 


where {//*} is defined as in Section 
the processes 


2.3 


Here, we furthermore define ig = Oi >4 ^ where 


v? = £ 


bi(t)dA 


(0 


>o Ci(t) + bi(t)A^l ( 
are assumed to be positive true ({7}, M)-martingales. 


Lemma 6.2. Let £ denote any non-negative 7 


( 1 , 2 , 3 ) 


-measurable random variable and let 


X = n ,>4 Xi where, for j = 4, 5, 6, Xj -measurable. Then 

Et R [tx]=K l^H^f [x7 

j> 4 

for M = M or Q and for 0 < t < T. 

Proof. Since 7^ 4,5,6 ^ is independent of 1:2 ' 3 ' 1 and of £, 


(6.46) 


E K 


d7 


( 1 , 2 , 3 ) 


V 7 


( 4 , 5 , 6 ) 


= e k 


d7 


( 1 , 2 , 3 ) 


Therefore, 

Ef Kx]=E R 
= E m 

= e m 
= e m 


e k 

e m 

e m 

t\F\ 


Cx|^ (1,2 ’ 3) V7^ 5 ’ 6) 


7 


( 1 , 2 , 3 ) 


v7 


( 4 , 5 , 6 ) 


£|7 

e|7; 

( 1 , 2 , 3 ) 


( 1 , 2 , 3 ) 


V 7 


( 4 , 5 , 6 ) 


x|7 


( 1 , 2 , 3 ) 


( 4 , 5 , 6 ) 


( 1 , 2 , 3 ) 


e k 


x 7 


( 1 , 2 , 3 ) 


V T\ 


( 4 , 5 , 6 ) 


x 17 


( 1 , 2 , 3 ) 


v 7 


( 4 , 5 , 6 ) 


v 7 


= Ef [e]Ef Lx] 


Next, the Girsanov formula in combination with the result for M-conditional expectation 
yields: 


«X']=Ef 

= Ef 


hTVrix 

IMVt 

VThTf, 

vtm 


= E: 


tt£, 


E? 


Tt 
TTffTX 

Tt”t 


E® 


VTX 

n 


= W K] e? Lx] • 


The result remains to be proven for the case £ = 1, which is done similarly. 
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For the XVA computations, we shall use a reduced-form counterparty risk approach in 


the spirit of Crepey ( 2012 ), where the default times of a bank “b” (we adopt its point of 


view) and of its counterparty “c” are modeled in terms of three Cox times r* defined by 


Ti = 


inf > 0 | jT 7 .W ds > E { 


(6.47) 


Under Q, the random variables E % (i = 4, 5, 6 ) are independent and exponentially dis¬ 
tributed. Furthermore, r c = 74 A 77 , = 75 A tq , hence r = 77 A r c = 74 A 75 A rg. 

We write 

c (4) . (6) b (5) . (6) (4) (5) . (6) 

It =Tt + Tt ,lt=Tt +Tt lit = 7t +Tt +Tt\ 

which are the so called ({J^}, Q)-hazard intensity processes of the {Gt} stopping times 
r c , Tb and r, where the full model filtration {Gt} is given as the market filtration {J 7 }- 


progressively enlarged by r c and Tb (see, e.g., Bielecki, Jeanblanc, and Rutkowski (2009), 
Chapter 5). Writing as before Dt = exp(— /(j r s ds), we note that Lemma 2.1 still holds in 
the present setup. That is, 

h t = ci(t) + bifCjA^ = D t fi t , 

an ({ J-) }, M)-supermartingale, assumed to be p ositive (e.g. under an exponential Levy mar¬ 
tingale specification for A^ ') as of Example 2.2). Further, we introduce = exp(— 7 ^ ds), 
for i = 4,5, 6 , and obtain analogously that 

k? ■= Ci(t) + bi(t)A^ = Z\ u[ l) . (6.48) 


With these observations at hand, the following results follow from Lemma 6.2 We write 


k t = Idi > 4 k {l) and z t = Ih>4 Z i 


W 


Proposition 6.1. The identities (2.22) and (2.26) still hold in the present setup, that is 

c 1 (T) + b 1 (T)A ( l 1) 


PtT = 

and, for t < Tj_ 1 , 


e - ft r ° ds 

= E® 

Dt 

= Ef 

hx 

L J 

[D t \ 


ht 


ci(t) + bi(t)A^ 


L(t\Ti-\,Ti) = 


L{ 0; T t - 1 , + b 2 {T i . 1 ,Tf)A^ ) + 63 ^- 1 , Tf)Af ] 


Likewise, 


E; 


e - ft T 7s ds 


= E; 


Zj* 

~Zf 


= E? 


Pot + b± (t)A t 

kx 


(i) 


h 


tt Cj(T) + bj(T)A ^ 
l=4,5,6 c i(t) T bi(t)Af 


(i) ’ 


(6.49) 


(6.50) 


(6.51) 


E; 


e -St 'Tsds^c 


= -El 


= -El 


Z, 


(5) 


L Z t 


(5) 


d T E) 


7(4) 7(6) 

Zj rj-i Zjrji 

L Z t (4)z ! 6 \ 


p -f t T 7 *dsl ST 


Ci(T) + bi(T)A\ 


(0 


i= 4,6 c *(^) + bi(T)A [ t 


(O’ 


(6.52) 
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(6.53) 


E; 


e ~ If (r a +^)ds 

— E? 

’ D t Z^zP' 


KJ t 

_ DtZ'fz f (6) _ 


tt Cj(T) + bj(T)A[ ^ 

i= 1 , 4,6 °i(t) + bi{t)A\ ■* 


Proof. Using Lemma |6.2[ we compute 

E? 


e ~ ft r ‘> ds 

= E? 

Dx 

= Ef 

JlT^T 

= Ef 

tlx 

Ef 

vx 

L J 


[D t \ 


. hm . 


. ht . 


. Pi _ 


= E; 


hx 


c i(T) + b\{T)A^ , 

ci(t) + b\{t)A^ 


(6.54) 


where the last equality holds by Lemma 2.1 This proves (6.49). The other identities are 
proven similarly. ■ 


Remark 6.3. Equations (6.49) and (6.51) are similar in nature and appearance. As it 


is the case for the resulting OIS rate {?’*} (2.17), the fact that (6.48) is designed to be a 


supermartingale has as a consequence that the associated intensity (6.45) is a non-negative 


process. This is readily seen by observing that } is a martingale and thus the drift 


of the supermartingale (6.48) is given by the necessarily non-negative process } that 


drives {zff. 

At time t = 0, all the Aq * = 0, hence only the terms Ci(T ) remain in these formulas. 


Since the formulas (6.49) and (6.50) are not affected by the inclusion of the credit component 


in this approach, the valuation of the basis swap of Section[5]remains unchanged. By making 
use of the so-called “Key Lemma” of credit risk, see for instance Bielecki, Jeanblanc, and 


Rutkowski (2009), the identity (6.53) is the main building block for the pre-default price 


process of a “clean” CDS on the counterparty (respectively the bank, substituting t/, for r c 
in this formula). In particular, the identities at t = 0 


3 “ Jo r ( r s+T'D' is 

fa (r s + Xs) ds 


= ci(T)c 4 (T)c 6 (T), 
= ci(T)c 5 (T)c 6 (T), 


(6.55) 

(6.56) 


for T > 0, can be applied to calibrate the functions cy(T), i = 4,5, 6, to CDS curves of the 
counterparty and the bank, once the dependence on the respective credit risk factors has 
been specified. The calibration of the “noisy” credit model components bi{T)A^\ i = 4, 5,6, 
would require CDS option data or views on CDS option volatilities. If the entire model is 
judged underdetermined, more parsimonious specifications may be obtained by removing 
the common default component (just letting r c = r 4 , Tf, = 75) and/or restricting oneself 
to deterministic default intensities by settting some of the stochastic terms equal to zero, 
i.e. bi(T)Af l ' > = 0, i = 4, 5 and/or 6 (as is the case for the one-factor interest rate models in 
Section [3]). The core building blocks of our multi-curve LIBOR model with counterparty- 
risk are the couterparty-risk kernels {k^}, i = 4, 5, 6, the OIS kernel {ht}, and the LIBOR 


kernel given by the numerator of the LIBOR process (2.26). We may view all kernels as 
defined under the M-measure, a priori. The respective kernels under the P-measure, e.g. 
the pricing kernel {7r^}, are obtained as explained at the end of Section [5] 
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6.2 XVA analysis 

In the above reduced-form counterparty-risk setup following Crepey (2012), given a contract 
(or portfolio of contracts) with “clean” price process {Pt} and a time horizon T, the total 
valuation adjustment (TVA) process {© 4 } accounting for counterparty risk and funding 
cost, can be modelled as a solution to an equation of the form 


©t = E' 


exp \ - J (r u + 7 U )du ) f s (@ s )ds 


, t £ [ 0 ,T], 


(6.57) 


for some coefficient {/)($)}■ We note that ( |6.57 ) is a backward stochastic differential 
equation (BSDE) for the TVA process {©t}- For accounts on BSDEs and their use in 
mathematical finance in general and counterparty risk in particular, we refer to, e.g., El 
Karoui, Peng, and Quenez (1997), Brigo et al. (2013) and Crepey (2012|) or (Crepey et al 


2014. Part III). An analysis in line with Crepey (2012) yields a coefficient of the BSDE 
(16.571) given, for d e M, by: 


m = 7 * c (1 - R c )(Pt - r*) + - 7 * (1 - R b )(Pt - r t )- 


CVA coefficient (cvat) DVA coefficient (dvat) 

+ b t T+ - b t Tf + A t (P t -d- r t ) + - A t{P t -d- r t )“, 

V -v-' 

LVA coefficient (lvat{d)) 


(6.58) 


where: 


- Rb and R c are the recovery rates of the bank towards the counterparty and vice versa. 

- r* = — Tj", where {T^} (resp. {r)~}) denotes the value process of the collateral 

posted by the counterparty to the bank (resp. by the bank to the counterparty), for 
instance Tj = 0 (used henceforth unless otherwise stated) or r t = P t . 

The processes {bt} and {bt} are the spreads with respect to the OIS short rate {rt} 
for the remuneration of the collateral {Tj 1 "} and {Tj7} posted by the counterparty and 
the bank to each other. 

The process {A*} (resp. {A*}) is the liquidity funding (resp. investment) spread of 
the bank with respect to {rt}- By liquidity funding spreads we mean that these are 
free of credit risk. In particular, 

A* = A t — (1 — Rb), (6.59) 

where {At} is the all-inclusive funding borrowing spread of the bank and where Rb 
stands for a recovery rate of the bank to its unsecured lender (which is assumed risk¬ 
free, for simplicity, so that in the case of {At} there is no credit risk involved in any 
case). 
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The data {T*}, {bt} and bt are specified in a credit support annex (CSA) contracted between 
the two parties. We note that 


E, 


exp ^ (r u + 7 U )duj f s (Q s )ds 


= E? 


= E? 


f WsDsZs 

It HsVtDtZt 


f s {Q s )ds 


h s k s „ \ , 

, , f s (®s)ds 

htkt 


(6.60) 


Hence, by setting 0 f = ht kt ©t, one obtains the following equivalent formulation of (6.57) 


and (6.58) under 


Q t = Ef 


^ fs(@s)ds 


(6.61) 


for t G [0, T] and where 


- A(j^) - * P - *■><* - r ‘> + - ^ - R ^ p ‘ - r ‘>' 


+ W-WT + M«-j^-r. 


- A ‘ |p ‘-ss- r ‘ 


(6.62) 


For the numerical implementations presented in the following section, unless stated other¬ 
wise, we set: 


7 b = 5%, 7 C = 7%, 7 = 10% 
Rb = R c = 40% 
b = b = X = X = 1.5%. 


(6.63) 


In the simulation grid one time-step corresponds to one month and m = 10 4 or 10' 5 scenarios 


are produced. We recall the comments made after (6.55) and note that (i) this is a case 
where default intensities are assumed deterministic, that is b{A 0) = 0 (i = 4,5,6) and 
(ii) the counterparty and the bank may default jointly which is reflected by the fact that 

7 1 < 7 1 + It- 


6.2.1 BSDE-based computations 


The BSDE (6.61)-(6.62) can be solved numerically by simulation/regression schemes similar 


to those used for the pricing of American-style options, see Crcpey, Gerboud, Grbac, and 


Ngor (2013), and Crepey, Grbac, Ngor and Skovmand (2014). Since in (6.63) we have 


A t = A t, the coefficients of the terms (Pt — — Ft) 41 coincide in (6.62). This is the case 


of a “linear TVA” where the coefficient ft depends linearly on d. The results emerging 


from the numerical BSDE scheme for (6.62) can thus be verified by a standard Monte 


Carlo computation. Table [l] displays the value of the TVA and its CVA, DVA and LVA 
components at time zero, where the components are obtained by substituting for d, in the 


respective term of (6.62), the TVA process ©t computed by simulation/regression in the 
first place (see Section 5.2 in Crepey et al. (2013) for the details of this procedure). The 
sum of the CVA, DVA and LVA, which in theory equals the TVA, is shown in the sixth 
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column. Therefore, columns two, six and seven yield three different estimates for ©o = 00- 
Table [2] displays the relative differences between these estimates, as well as the Monte Carlo 
confidence interval in a comparable scale, which is shown in the last column. The TVA 
repriced by the sum of its components is more accurate than the regressed TVA. This 
observation is consistent with the better performance of Longstaff and Schwartz (2001) 
when compared with Tsitsiklis and Van Roy (2001) in the case of American-style option 
pricing by Monte Carlo methods (see, e.g., Chapter 10 in Crepey (2013)). 


m 

Regr TVA 

CVA 

DVA 

LVA 

Sum 

MC TVA 

10 4 

0.0447 

0.0614 

-0.0243 

0.0067 

0.0438 

0.0438 

10 5 

0.0443 

0.0602 

-0.0234 

0.0067 

0.0435 

0.0435 


Table 1: TVA at time zero and its decomposition (all quoted in EUR) computed by re¬ 


gression for m = 10 4 or lO 5 against and X^'. Column 2\ TVA ©o- Columns 3 to 5: 
CVA, DVA, LVA at time zero repriced individually by plugging for $ in the respective 
term of (6.62). Column 6: Sum of the three components. Column 7: TVA computed by a 
standard Monte Carlo scheme. 


( 2 ) 


m 

Sum/TVA 

TVA/MC 

Sum/MC 

Cl// MC 

10 4 

-2.0114% 

2.0637% 

0.0108 % 

9.7471% 

10 5 

-1.7344 % 

1.7386 % 

-0.0259% 

2.9380% 


Table 2: Relative errors of the TVA at time zero corresponding to the results of Table [I] 
“A/B” represents the relative difference {A — B)/B. “CI//|MC|”, in the last column, refers 
to the half-size of the 95%-Monte Carlo confidence interval divided by the absolute value 
of the standard Monte Carlo estimate of the TVA at time zero. 


In Table [3j in order to compare alternative CSA specifications, we repeat the above 
numerical implementation in each of the following four cases, with A t set equal to the 


constant 4.5% everywhere and all other parameters as in (6.63): 


1. 

( Rb, Rb, Rc) 

= (100,40,40)%, 

Q = P, 

r = 

2. 

(Rb, Rb, Rc) 

= (100,40,40)%, 

Q = P, 

r = 

3. 

(Rb, Rb, R c ) 

= (40,40,40)%, 

Q = P, 

r = 

4. 

(Rb, Rb, R c ) 

= (100,100,40)%, 

Q = P, 

r = 


(6.64) 


Remembering that the t = 0 value of both legs of the basis swap is equal to EUR 27.96, 
the number in Table [3] may seem quite small, but one must also bear in mind that the toy 


model that is used here doesn’t account for any wrong-way risk effect (see Crepey and Song 


(2015j)). In fact, the most informative conclusion of the table is the impact of the choice of 
the parameters on the relative weight of the different XVA components. 
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Case 

Regr TVA 

CVA 

DVA 

LVA 

Sum 

Sum/TVA 

1 

0.0776 

0.0602 

-0.0234 

0.0408 

0.0776 

-0.0464 % 

2 

0.0095 

0.0000 

0.0000 

0.0092 

0.0092 

-3.6499% 

3 

0.0443 

0.0602 

-0.0234 

0.0067 

0.0435 

-1.7344 % 

4 

0.0964 

0.0602 

0.0000 

0.0376 

0.0978 

1.4472% 


Table 3: TVA at time zero and its decomposition (all quoted in EUR) computed by regres¬ 


sion for m = 10 5 against X\ and X\ . Column 2: TVA ©o- Columns 3 to 5: CVA, DVA 
and LVA at time zero, repriced individually by plugging @t for d in the respective term of 
(6.61). Column 6: Sum of the three components. Column 7: Relative difference between 
the second and the sixth columns. 


6.2.2 Exposure-based computations 

Let’s restrict attention to the case of interest rate derivatives with {Pt} adapted with respect 
(1 2 3 ) 

to {J-f ' ’ }■ We introduce c(s) = rii> 4 c *(s) an d the function of time 

EPE(s) := E M [h s P+] = E Q [D s P+] , resp. E M [h a P~] = E Q [D 9 P~] , 


called the expected positive exposure, resp. expected negative exposure. For an interest- 
rate swap, the EPE and ENE correspond to the mark-to-market of swaptions with maturity 
s written on the swap, which can be recovered analytically if available in a suitable model 
specification. In general, the EPE/ENE can be retrieved numerically by simulating the 
exposure. 

In view of (6.61)-(6.62), by the time t = 0 forms of (6.51) and (6.52), the noncollater- 
alised CVA at t = 0 satisfies (for R c ^ 1, otherwise CV Aq=0)\ 


(1 -Rc 


-CVA 0 = e r 


h s k s ^ c s P+ ds 


= / E R 


uo 
[h s P+] 1 


= / E m [h 8 P+] E m [k s Y s ] ds 
c 6 (s 


'0 


[Z s7s c ] ds = - EPE(s 

Jo 


+ 


C 4 (S, 


k c 6 (s) c 4 (s) 

Similarly, for the DVA (for Rf, ^ 1, otherwise DVAq = 0) we have: 

1 


c(s)ds. 


(l - R b ) 


DVAq = - 


ENE(s) 


L^ + mM |c(s) *. 

C_ 2 (s) C_l(s) 


For the basis swap of Section [5] and the counterparty risk data (6.63), we obtain by this 
manner CV Aq = 0.0600 and DVAq = —0.0234, quite consistent with the corresponding 
entries of the second row (i.e. for rn = 10 5 ) in Table [lj As for the LVA, to simplify its 
computation, one may be tempted to neglect the nonlinearity that is inherent to Ivat(’d) 
(unless Xt = A t), replacing i? by 0 in lvat(’&). Then, assuming lvat( 0) e , by (6.57)- 

(6.58), one can compute a linearised LVA at time zero given by 

rT 


LVA 0 = E K 


h s k s lva s (0)ds 


LJ0 


_ j ]g M \j lg ^a s (0)] E m [A: s ] ds = f E M [h s Zi;a s (0)] c(s)ds, 


io 
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by (6.51) for t = 0. This is based on the expected (linearised) liquidity exposure 


E M [h s lva s { 0)] = E Q [£>>a s (0)] . 

In case of no collateralisation (Tt = 0) and of deterministic A t and A t, we have 

lva s {0) = X S P+- X S P~, LVA 0 = [ (x s EPE{s) - X s ENE(s)) c(s)ds. 


In case of continuous collateralisation (Tf = Pt) and of deterministic bt and bt, the formulas 
read 

lva s (0) = b s P+ -b a P~, LVA 0 = [ (b s EPE(s) — b s ENE(s)) c(s)ds. 

Jo 

As for CVA/DVA, the LVA exposure is controlled by the EPE/ENE functions, but for 
different “weighting functions”, depending on the CSA. For instance, for the data (6.63), 
the LVA on the basis swap of Section [ 5 ] (collateralised or not, since in this case bt = h = 
Af = Af = 1.5%), we obtain LVAq = 0.0098, quite different in relative terms (but these are 
small numbers) from the exact (as opposed to linearised) value of 0.0067 in Table fll We 
note that the stochasticity of the default intensities ( 7 )* } is averaged out in all these t = 0 
pricing formulas, but it would appear in more general t pricing formulas or in the XVA 
Greeks even for t = 0. 
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